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Abstract 

It is shown that the relativistic jet, emitted from the center of the Galaxy 
during its activity, possessed power and energy spectrum of accelerated protons 
sufficient to explain the current cosmic rays distribiition in the Galaxy. Proton 
acceleration takes place on the light cylinder surface formed by the rotation 
of a massive black hole earring into rotation the radial magnetic field and the 
magnetosphere. Observed in gamma, x-ray and radio bands bubbles above and 
below the galactic plane can be remnants of this bipolar get. The size of the 
bubble defines the time of the jet's start, ~ 2.4 • 10'' years ago. The jet worked 
more than lO'' years, but less than 2.4 • 10^ years. 
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1. Introduction 

The traditional point of view on the; origin of cosmic rays in the Galaxy is the 
concept of acceleration of charged particles at fronts of shocks from supernova 
explosions. Arguments in favour of this mechanism are sufficient mechanical 
energy that is released when the supernova explodes, as well as universal index 
of power law spectrum of particles, accelerated by strong shocks. Total power 
of cosmic ray sources in order to maintain their observed density of energy is 
5- 10**^ er(//s, which equals approximately 15% of the kinetic energy of supernova 
explosions. When the gas compression in a shock is equal to 4, the index of the 
power law energy spectrum of accelerated particles is equal to -2, N{E) oc 
It is in a good agreement with observed cosmic ray spectrum at energies 
£ < 3 • W^eV. Beginning from the first works by Krymskii, 1977; Bell, 1978; 
Blandford & Ostriker, 1978, who proposed the mechanism of acceleration of 
charged particles on fronts of shocks propagating in the turbulent environment, 
much progress has been made to explain the observed characteristics of galactic 
cosmic rays in the belief that they are accelerated at the front of shocks. 

On the other hand, there is no objections to generate galactic cosmic rays 
in one source in the Galaxy (Ptuskin & Khazan, 1981). This potential source 
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can be the center of the Galaxy, which is the massive black hole of M ~ 4 • 
IO^Mq mass. And while the luminosity of Sgr A* is small now, it is only 10^^ 
erg/s, in the past the center could be much brighter because its Eddington 
luminosity equals LEdd = 5.2 • 10'*^ crg/s. On the past activity of the center of 
the Galaxy shows newly discovered by Fermi Gamma-ray Space Telescope above 
and below the galactic plane big bubbles emitting gamma radiation in the range 
of 0.1 — 1000 GeV (Su et al., 2010). Such formations was previously observed 
in the x-ray range (1.5 - 2) KeV by ROSAT All-Sky Survey (Snowden et al, 
1997) and in the microwave range (20-40) GHz by WMAP (Finkbeiner, 2004). 
Estimated energy stored in bubbles is of 10'"'^ — 10^^ erg (Sofuc, 2000). As we 
will see below, bubbles of a relativistic gas could be formed by the jet, emitted 
from surroundings of the massive black hole. Here we provide an alternative 
mechanism of origin galactic cosmic rays, in which the nucleus of the Galaxy 
in the active phase injected the relativistic jet, which was the source of cosmic 
rays. 

In the following sections we will calculate the power of the jet and the energy 
spectrum of protons in the relativistic jet, as well as describe the remnants of 
the relativistic jet injected from the center of the Galaxy, having the form of 
bubbles above and below the galactic plane. In the final section we will discuss 
correspondence of the scenario of the galactic cosmic rays origin provided with 
cosmic rays characteristics observed. 

2. Relativistic jet 

Sources of energy of active galactic nuclei are the accretion on a massive black 
hole, in which the gravitational energy of a falling gas transforms into radiation 
and heat, as well as the rotation of a black hole. Mechanism of extraction of 
energy and angular momentum from the black hole is called as the mechanism of 
Blandford-Znajek (1977). The energy of a rotating black hole is a large value, 
Erot = = a^Mc?/8 = 2.25 • 10^^a^{M/MQ) erg. For the Galaxy 

Erot — 9 • 10'''^a^ erg. Here we have introduced the dimcnsionless parameter of 
a, describing rotation of the black hole, a = Jc/]VPG, a < 1. J is the black 
hole angular momentum, G is the gravitational constant. Angular velocity of 
rotation of a black hole is proportional to the value of a, = ac/'^rH, Th 
is the gravitational radius of not rotating black hole (a = 0), r// = 2MG/c?. 
Energy extraction is possible when there is a poloidal magnetic field B near the 
black hole horizon. In this case, rotating black hole acts as a Dynamo machine, 
creating a voltage U = fn^H/'^T^c (Landau & Lifshits, 1984). The value of 
fn is the flux of the poloidal magnetic field reaching the horizon of a black 
hole, Jh — TrBrfj. Voltage U generates the electric current / = U/{R + Rh), 
which on the one hand is closed on the horizon of a black hole that has the 
resistance Rh = 47r/c ~ 377 ohm. Resistance of the outer part of the current 
loop is R. Thus, the power, extracting from a rotating black hole, is L = RP = 
U'^R/{R + RhY = a^B^rjjR/16{R+ Rh)^, and reaches the maximum L^, at 
R = Rh, Lm — B^r'^c/ 25671. The value of L is proportional to the energy 
of the poloidal magnetic field near a black hole and can reach the Eddington 
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Figure 1: Configuration of the magnetic field and electric currents in the jet and in the disk. 

luminosity at sufficiently large magnetic fields B ~ 10^a~^ Gauss in the center 
of the Galaxy. This field is accumulated near the horizon of a black hole in the 
process of accretion of a disk matter in which the magnetic field is frozen. Thus, 
for the effective work of the mechanism of Blandford-Znajek an accretion disk 
around a massive black hole is required, not as a source of the energy, but as 
the agent bearing the magnetic field to a black hole. In addition, the electric 
current / fiows in the disk, this is the part of the current loop: in the disk, in the 
black hole horizon, then in the jet, closing at large distances in the interstellar 
matter (see Figure 1). 

In the disk, in addition to the radial electric current Ip = J jpds = —I, 
stronger toroidal current j^, ~ lO^jp, flows also (Istomin & Sol, 2011), it 
generates the radial magnetic field B. 

The rotating black hole brings into rotation the radial magnetic field in the 
magnetosphere of a black hole above and below the disk. Angular velocity of 
rotation of the magnetic field lines, the same as rotation of the magnetospheric 
plasma, fJi? , is proportional to the angular velocity of rotation of the black hole 
Hh, = ^hR/{R + Rh) (Thorne et al., 1986). Plasma rotation is the drift 
motion in crossed radial magnetic field and electric field of plasma polarization. 
Thus, there appears so-called the light cylinder surface in the black hole magne- 
tosphere, where the magnitude of the electric field is compared with that of the 
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magnetic field and the rotation velocity approaches the speed of light c. The 
radius of the light surface is vl = c/flp = 2a~^rH{R + Rh)/R > th- On the 
light surface charged particles get considerable energy and angular momentum 
of rotation. Energy density of particles on the light surface is compared with the 
energy density of the electromagnetic field (£'£ + B\)/Sit = Bj^/An (Istomin, 
2010). All energy passes to protons, 7 = Sp/nipC^ >> 1. Electrons are practi- 
cally not accelerated due to large synchrotron losses in a strong magnetic field 
(Istomin & Sol, 2009). Energetic protons, accelerated near the light surface, 
and whose energy is mainly in the azimuthal motion, create the jet. Jet's power 
is Lj = B'^r'j^c{u}cHrH /c)~^^'^ /2 (Istomin & Sol, 2011). Here cVcH is the non 
relativistic cyclotron frequency of protons in the magnetic field near the black 
hole, cOcH = eB/nipC. Jet arises when the power extracted from the rotating 
black hole L becomes greater than the jet power Lj, L > Lj. This imposes a 
limitation on the value of the magnetic field 



> (1287r)^a-« 



[R + Rh? 



4:RRh 



(1) 



For R = Rh it gives 



B > 2.7- 10"a-«^Gauss. (2) 

M 

For the center of the Galaxy, the magnetic field must satisfy the condition 
B > 6.75 • 10^a~^ Gauss. We see that to generate a jet, less massive black holes 
should have a stronger poloidal magnetic field near the horizon, B oc M~^. In 
addition, rotation must be fast, close to the critical value of a ~ 1, because of 
strong dependence of the expression (2) on a. It should also be noted that the 
resistance of the external current loop i?, on which the jet power Lj depends, 
is not the ohmic one Rc, which turns out to be small, Rc « Rh (Istomin & 
Sol, 2011), but is the effective resistance Rj, which can be attributed to the jet, 
receiving energy from the rotating black hole. If L = Lj the resistance of the 
jet is Rj = Rh- Under the equality in the expression (1) when a rotating black 
hole can begin to generate a jet, the expression for the jet's power becomes 
universal 



Lj = 2*«7r^mpc^ ) "~ = ^'^ ' 10 ^r^/s. (3) 

All jet energy are in the energy of protons, and, as we can see, is sufficient for 
the production of galactic cosmic rays. 



3. Energy spectrum of fast peirticles in jet 

Istomin and Sol (2009) had shown that on the light surface, produced by 
the rotating radial magnetic field, which is carried into rotation by a black 
hole, protons gain considerable energy. The Lorentz factor 7 becomes equal 
to 7 = (7o7i)^^^- The value of 7j is the Lorentz factor of particles in the 
magnetosphere of a black hole before crossing the light surface. And the value of 
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7o is the maximum of the Lorentz factor, which could be achieved by a particle in 
this acceleration mechanism, 70 = Wcl/^f- Here oJcL is the cyclotron frequency 
of rotation of protons in the poloidal magnetic field near the light surface. When 
7 = 70 the cyclotron radius of a proton is compared with the radius of the light 
surface. For non relativistic particles of the black hole magnetosphere, 74 ~ 1, 
the Lorentz factor of accelerated particles is equal to 7 = 7^^ = (wcl/^f)^^^- 
Crossing the light surface at different distances z from the accretion disk plane, 
particles gain different energies, since the magnetic field decreases with distance 
from the black hole. For a radial magnetic field B oc (z^ + r\)^^. Thus, 
7 — 7m(l + /r\)~^/'^ , where 7^ is the maximal Lorentz factor of accelerated 
particles near the accretion disk. Accelerated protons of the jet are collected 
from various parts of the light cylinder surface of rj, radius, but located at 
different distances z. Therefore, the number of particles is dN oc ndz, where n 
is the density of protons in the magnetosphere near the light surface. Connecting 
values z and 7, we get 

dz_ 7m (i7 

Considering that the vertical size of the magnetosphere is larger than the light 

surface radius, the density n can be taken as constant. As a result we get the 
distribution function of relativistic protons in the jet, ^(7) = dN/d'y, 

F(7) = const • 7-'(l - 7 < 7m- (4) 

Wc see that in the range 7 << 7™ the spectrum of relativistic protons is the 
power law spectrum with the index -2. This spectrum is observed in gamma 
radiation from bubbles above and below the Galactic plane by Fermi Gamma- 
ray Space Telescope (Su et al., 2010). Considering that the gamma radiation 
occurs due to collisions of relativistic protons with interstellar gas through meson 
production (Crocker & Aharonian, 2011), and the distribution of photons is 
similar to the distribution of protons, one can conclude that the jet from the 
center of the Galaxy actually existed, and bubbles are filled with relativistic 
protons of the jet. The value of 7^ is (Istomin & Sol, 2011) 

7m=(^j ' 

and for R = Rh equals (see the expression (1)) 

7™ = [USnfa-^ = 1.6 • lO^a"^. 

The spectrum of protons (4) breaks at 7 = 7™ and has there the root singularity 
(integrable) that is smoothed considering the thermal dispersion of particles in 
the magnetosphere of the black hole, Aji = Sp/nipC^. The distribution (4) is 
shown on Figure 2. The value of 7m is chosen to be equal to the Lorentz factor 
of the break in the observed spectrum of cosmic rays at the energy = 3 • 10^'' 
GeV, 7to = 3.2 • 10^. This corresponds to the rotation parameter a = 0.47. The 
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Figure 2: Distribution function of relativistic protons in the jet, 7 < 7m- The value of 7, 
corresponds to the brealc in the spectrum of the cosmic ray. The slope is equal to -2. 



power of the jet (3) is Lj ~ 8.9 • 10^^ erg/s. That is in agreement with estimated 
from observations powers of jets ejected from active galactic nuclei-10'''^ — 10'''^ 
erg/s (Mao-Li et al., 2008). 

Relativistic protons with the spectrum (4) are formed from thermal par- 
ticles of the black hole magnetosphere, 7^ ~ 1. But in addition to thermal 
particles in the magnetosphere there can exist accelerated protons. The tur- 
bulent motion of the accreting disk matter in the presence of the frozen mag- 
netic field leads to acceleration of particles, which have the power law energy 
spectrum, f{-f) = const ■ 7"*^, 7 < 71, /3 ~ 1 (Istomin & Sol, 2009). The 
disk must bo turbulent to provide for the abnormal gas transport. Getting 
onto the light surface, accelerated protons are converted to more energetic, 
7 — > (707)^^^- Their distribution function becomes equal (Istomin & Sol 2009) 

f'{j) = 2const ■ 7^"V"^^"'^S 70^^ < 7 < (7071)^/^- Thus, there is another 
component of jet relativistic protons, their number is equal to 

iVoc/ dzl^^ 7o''"V^+'rf7,7o=7^(l + ^>i)-^ (6) 



Jjg' 



Transforming the integration area in (6), we get 

N oc / ^-^f^+^d-/ / 1 + ^ — + 



The first term in Eq. (7) corresponds to relativistic protons with energies 7 < 7to 
similar to protons (4), accelerated from the thermal gas. Their distribution 
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Figure 3: The distribution function of relativistic protons in the jet, 7 > 7m- The value of 
7m corresponds to the break in the cosmic ray spectrum. The spectral index equals -2.4. The 
maximum energy is 10^® eV. 



function equals 



F(7) = coast ■ 7-2^+1 / 



In the energy range 7 << j„i this distribution has the same power law spec- 
trum (4), F{j) oc 7^^. But since the number of accelerated particles in the 
magnetosphere of the black hole is much less than that of thermal particles, 
the contribution of these particles into the total distribution at 7 < jm can be 

neglected. The second term in Eq. (7) describes the distribution of relativistic 

1/2 

protons at 7 < 7to < 7m7i 



(75.7i/7'-l)'^' / .2xl-f* 



/■(7m 71/7 -1) ' / 

F(7) = const ■ 7-2^+1 J (^1 + 



' L 



dz 



\const ■ 7-^''+^ [B{1, /3 - 3/2, 1/2) - i3(7V(7^7i), P " 3/2, 1/2)] . (8) 

Here B{x^ a, b) = t'^~^{l — t)''~^dt is the incomplete Beta function, B{1, (3 — 
3/2, 1/2) = Tr^/^r{P - 3/2)/r(/3 - 1), r(a;) is the Gamma function. The distri- 
bution (8) with ^ = 1.7, 7m = 3.2 • 10^ and 71 = 10^ is shown on Figure 3. For 

1/2 

energies 7 < 7m 7i the distribution of relativistic protons is the power law with 
index —(2/3 — 1). At 7 ~ 7m7y^ the distribution falls down, for 7^, = 3.2 ■ 10^ 
and 7i = 10^ the maximum energy is ~ 10^* eV. 

We have chosen the value of /3 = 1.7 from the fact that indices of spectrum 
of cosmic rays before and after the break at energy 3 • lO'^^ eV differ on the 
value of 0.4 - the spectrum becomes more soft with the index —3.1. The same 
difference in indices must be in the source of cosmic rays, which in our case 
is the relativistic jet. At /3 = 1.7 the index of the spectrum (8) equals -2.4, 
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while the index of the speetrum (4) is equal to -2. It should be noted that the 
distribution function of relativistic protons (4) at energies 7 < 7^, and (8) at 
energies 7 > 7™, is continuous, i.e. F{j = 7^ — 0) = ^(7 = 7m + 0). This is 
because the acceleration on the light surface undergo protons of the black hole 
magnetosphere with a unique spectrum - thermal at low energies, turning into 
the tail of fast particles up to the energy 7 = 71. If their energy distribution 
function is Fi{£i), then the distribution of particles, accelerated on the light 
surface, have also the continuous distribution F{£) = Fi[£i{£)]d£i/d£,£i = 
£'^/^omp(?. 

Here and hereafter, we talk about protons, bearing in mind that they are 
'heavy' particles, unlike electrons. Nuclei of m mass and Ze charge will be 
accelerated effectively also on the light surface. As we saw above, the efficiency 
of acceleration depends on the value of 7^^ = [lOc/^fY^'^ ■ Thus, nuclei will 
receive energy per nucleon [Zmp/mY^'^ ~ 2~^/^ times less than protons. 

4. Jet's remnants 

The center of the Galaxy, being active, created a relativistic jet, particles 
of which spread in the Galaxy, forming isotropic background of cosmic rays. If 
one consider that the angular momentum of the massive black hole in the ccmter 
of the Galaxy coincides with that of the Galaxy, than the direction of the jet 
propagation is perpendicular to the plane of the Galaxy. The jet length can be 
quite large, so jet in M87 extends ~ 2 kpc. Therefore, above and below the 
Galactic plane (if we have two almost symmetrical jets) one can see traces of 
the jet. 

Let us consider a simple diffusive model of propagation of relativistic parti- 
cles, whose source is located in the center of the Galaxy (r = 0), in an interstellar 
medium 

^ + uVTV - = Q{t)5{T). (9) 

Here u is the velocity of the interstellar matter. We consider the region outside 
the stellar galactic disk, then the velocity u is the speed of the galactic wind, 
which is along the coordinates z which is perpendicular to the plane of the 
Galaxy. The; value of D is the diffusion coefficient, which can be anisotropic. 
Diffusion depends on the intensity of the magnetic field B, falling exponentially 
with the distance z from the galactic plane, B = Bq exp(— z/zi), zi ~ 2 kpc. 
Diffusion of charged particles is determined by their motion in the magnetic 
field, and the diffusion coefficient is inversely proportional to the magnetic field 
strength, D oc B~°'. So, for the most strong Bohm diffusion D = CTcji-, Tc 
is the proton cyclotron radiiis, a = \. Thus, the particle diffusion increases 
exponentially with the coordinate z, D = exp(z/2;o), zq = Zx/a. Such a 
strong dependence of the diffusion on the coordinate z leads to effective non 
diffusion expansion of particles along this coordinate and decreasing its density. 
To take into account this effect we insert the function (p = Nex.p{z/zo). Eq. 
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(9) becomes 



exp ( f + «cxp (-^) f ^ - ^) + 

' zoj dt \ zqJ \dz zo J zo dz 



Dn 



p dp \^ dp J ' dz 



K d / d(p\ d'^ifi'^ 



We consider the distribution of particles as azimuthal symmetric depending 
on the distance z and the cyhndrical radius p. The value of k is the ratio of 
the transverse diffusion coefficient D±, perpendicular to z, to the longitudinal 
one Dy, along z, k = Dj_/D\\. We see that in the equation of particle motion 
there appears the effective velocity along z, uq — Do/zo. It arises as a result 
of the exponential growth of the particle diffusion over z. For characteristic 
values of -Do = 5 • 10^* cm? /s and zq = 2 kpc the velocity uq is of uq — 10^ 
km/s. This speed is much larger than the speed of the galactic wind u, which at 
distances of several kpc from the galactic plane is less than 30 km/s (Ptuskin, 
2007). Although the wind velocity increases with distance z, the exponential 
factor exp(— z/zq) in Eq. (10), allows us to ignore the velocity of the galactic 
wind in comparison with the velocity mq. Values of Dq and zo specify scales of 
length and time, so it is convenient to go to the dimensionless variables in Eq. 
(10) z' — z/ Zq, p' = p/zo, t' = t{Do/zQ), ip' = (fZQ. Eq. (10) becomes (primes 
are omitted) 

exp(-.)$ + ^ - ft - (p^) = Qit)'-^. (11) 



dt dz dz"^ p dp \ dp J 2irp 

At z > 1 the propagation (first derivative over z) predominates over the diffusion 
(second derivative), and Eq. (11) makes easy 

^ ^ Q[t-{l-e-^)] 

Attkz \ Akz 

From this solution one can see that to the point z come particles, which were 
emitted by the source at the retarded time t' = t — {1 — e^). If the time of 
the jet's start is t = 0, particles will lift to the maximum height of Zm = 
— ln(l — t), < < 1. Formally at t = 1 particles come to infinity during the 
finite time, that is impossible. The velocity limitation implies the condition 
Zm < {czo/Do)t, c is the speed of the light, which is not difficult to hold because 
czq/ Dq ~ 3.6 • 10'^. Knowing the value of Zm from observation one can estimate 
the time t\ when the power source of cosmic rays in the center of the Galaxy 
begins to work, i.e. when the jet starts, = 1 — exp(zm). In dimensional units 
ti = to[l - exp(zm/zo)], to = Zq/Dq. When Zq = 2 kpc and I?o = 5 • 10^^ err? /s 
the time is Iq = 7.6 • 10^* s = 2.4 • 10^ yr. The gamma radiation observed 
above and below the galactic plane extends to the height of about 8 kpc (Su 
et al., 2010), i.e. Zm — 4. This means that in fact ti = to. In addition, 
one can find the time t2 when the source turns off. If the jet worked a short 
time, all particles would have lifted in height at z > 1, and we would see their 
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Figure 4: The density distribution of relativistic particles (12) above the galactic plane in 
dimensionless coordinates: 2, the distance from the plane of the Galaxy, and p, the cylindrical 
radius. 



absence near the galactic plane. Since this is not observed in bubbles, then 
t2 < toil - 1/e) ~ 0.6^0- 

Knowing the solution for 1^9 at z > 1 we can find the density of relativistic 
particles N{t, z, p) in the same region 



Q[t - (1 - e--)] _ \ f 



N = — — g exp 



h z + ln(z' 



(12) 



The distribution N{z,p) (12) for the permanent source, Q — const{t), is shown 
on Figure 4. We also draw levels of the constant density N, p^ / Akz + z + ln{z) = 
const, Figure 5. The value of k, the anisotropy of diffusion, is chosen to be equal 
to K = 0.14 from those considerations that the observed ratio of bubble's scales 
z/ p ~ 8kpc/3kpc = 8/3 would be consistent with the forms of the constant 
density profiles painted on Figure 5. It seems that such value of k indicates that 
the magnetic field in the halo of the Galaxy near the center is mostly vertical. 
And this is natural because for the cylindrical symmetry radial and azimuthal 
magnetic field components should approach zero on the axis p = 0. 



5. Discussion 

We have shown that whereas in the past the center of the Galaxy was active 
and radiated the jet, its energy and composition are sufficient to explain the 
origin of cosmic rays in the Galaxy. Bubbles of a relativistic gas observed in 
gamma, x-ray and radio bands above and below the galactic plane, apparently, 
are remnants of the bipolar jet emitted from the vicinity of the massive black 
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Figure 5: Levels of the constant density of relativistic particles specified by the distribution 
(12). 
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hole in the center of the Galaxy. The vertical size of the bubble z ~ 8 kpc 
permits us to estimate the time of switching on of the jet, ~ to = 2.4- 10 years 
ago. We can also estimate the lower limit of the jet's work, At = t\ —t2 > OAto — 
10^ yr. During the time At the jet got the energy LjAt ~ S.Q-IO"*^ erg/s x 3-10"'^'* 
s ~ 2.7 • 10^° erg, that is slightly larger than the energy stored in the black hole 
rotation ~ 10^*^ erg. However, if we estimate the mass of the accreted matter, 
absorbed by the black hole at the same time, it can reach a large part of the 
mass of the black hole, AM ~ MEdd^t = 9.2 • IQ-^Mq yr''^ x 10^ yr ~ IO^Mq, 
AM ~ M/4. Transmitted to the black hole by the accreted matter, the angular 
momentum AJ can even exceed its initial value (Istomin, 2004). The jet's 
energy ~ 10^° erg is enough to fill by cosmic rays as the disk (10^^ erg), as 
the halo (10^'^ - 10^^ erg) of the Galaxy. Filling of the disk of the Galaxy by 
relativistic particles is described by the same Eq. (11), but there should be 
l^l < 1. Therefore, we can ignore the dependence of the diffusion coefficient on 
the distances z, and solve the pure diffusion equation 

- Dw^Tn- - D±-^^ = Q{t)- 



dt " dz^ pdp\ dp J 2-jTp 

Here we are interested in distribution of particles in the disk on the transverse 
distances p, so we average this equation over z and get 

where the value of D± = Dg is the diffusion coefficient of cosmic rays in the 
galactic disk. The solution is 

The time ti is the start time of the jet. If the jet had worked with constant 
power Q and switched off at the time t^, then the density distribution of cosmic 
rays in the disk is 

N 



/ x-^e-'^dx. (14) 



The solution (14) shows that if jADgti < 1 then the distribution of the density 
over the radius p is almost uniform. For p^ /ADg{ti — t2) > 1 it is logarithmic, 
N (X - \n{p^ I'iDgti), and for p^ / ADg{ti-t2) < 1 it is constant, N oc \n{ti/{ti - 
^2)). The diffusion coefficient of cosmic rays in the disk equals Dg = 2.2 • 
^028^0.6 g^2/g (ptuskin, 2007). The condition R'^/ADg <to, Rc^ 15 kpc is the 
disk radius, imposes a lower limit on the energy of protons, homogeneously filling 
the galactic disc, 7 > 400. Apart the diffusion particles can move freely along 
the regular magnetic field of spiral arms of the galactic disk. The necessary 
for filling velocity R/to — 6 • 10 cm/s = 2 • 10"'^ c does not contradict the 
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observed anisotropy of cosmic rays S, S ^ lO"''. Since the dependence of the 
distribution (14) over the energy is determined not only by the energy spectrum 
of the source Q(7) but also the dependence of the diffusion coefficient Dg oc 7°'^ 
over the energy, the spectrum of particles in the disk will be softer than that 
in the source. Thus, the discussed mechanism of origin of galactic cosmic rays 
by the jet, emitted from the center of the Galaxy, satisfactorily explains the 
observed spectrum, the index -2.7 (-2.6 for the jet) before the break, and the 
index -3.1 (-3.0 for the jet) after the break. 

Cosmic rays, filling simultaneously the galactic disk and the halo, fiow out 
from the Galaxy. Their lifetime r is determined by as the energy loss time te, 
as the time of diffusion leakage from the Galaxy after the source of relativistic 
particles stopped, Tn- The time te is estimated as — 3 • 10^ yr (Strong & 
Moskalenko, 1998). It is larger than the time of the jet's start ti, i.e. beginning 
of the filling of the Galaxy by cosmic rays, te > ti. The diffusion time td = 
/AD is determined by the distance r that particles travel during the period 
from the switching on of the source r = (4Dti)-'^/^. Thus, td ~ ti does not 
depend on the energy of particles. This time is also larger than the time of jet's 
switching off t2, i2 < 0.6fi. We see that to now the distribution of relativistic 
particles, generated by the jet, does not change noticeable. 
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